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We present a model for compressible sub-Alfvenic isothermal magnetohydrodynamic (MHD) tur- 
bulence in low j3 plasmas and numerically test it. We separate MHD fluctuations into 3 distinct 
families - Alfven, slow, and fast modes. We find that, production of slow and fast modes by Alfvenic 
turbulence is suppressed. As a result, Alfven modes in compressible regime exhibit scalings and 
anisotropy similar to those in incompressible regime. Slow modes passively mimic Alfven modes. 
However, fast modes show isotropy and a scaling similar to acoustic turbulence. 

PACS numbers: 52.35.Bj, 47.65. +a, 52.30.-q 52.35.Ra 



I. INTRODUCTION 



Most astrophysical fluids, including stellar winds and 
the interstellar medium (ISM), are turbulent El |7j] with 
an embedded magnetic field that influences almost all 
of their properties. High interstellar Reynolds numbers 
(Re = LSV/v > 10 8 ; L=the characteristic scale or driv- 
ing scale of the system, <5V^=the velocity difference over 
this scale, and v— viscosity) ensure that. Turbulence 
spans from km to kpc scales and holds the key to many 
astrophysical processes (e.g., star formation, fragmenta- 
tion of molecular clouds, heat and cosmic ray transport, 
magnetic reconnection). Statistics of turbulence is also 
essential for the CMB foreground studies jlf| . 

Kolmogorov scalings |D| were the first major advance 
in the theory of incompressible (non-magnetized) turbu- 
lence. Kolmogorov theory predicts an isotropic power law 
energy spectrum (E(k) oc /c~ 5 / 3 ) in wave- vector space k. 

Attempts to describe magnetic incompressible turbu- 
lence statistics were made by Iroshnikov [ ll| and Kraich- 
nan |l4j. Their model of turbulence (IK theory) is 
isotropic in spite of the presence of the magnetic field 
and predicts fc -3 / 2 power law energy spectra for both 
velocity and magnetic field. However, the assumption 
of isotropic energy distribution in wave-vector space has 
been criticized by many researchers p2| . 

An ingenious model very similar in its beauty and sim- 
plicity to the Kolmogorov model has been proposed by 
Goldreich & Sridhar || (hereinafter GS95) for incom- 
pressible MHD turbulence. It predicts a Kolmogorov-like 
energy spectra (E(k±) oc k ± 5 ^ 3 ) in terms of wave- vector 
component which is perpendicular to the local direc- 
tion of magnetic field. The parallel component of the 

2/3 • • 

wave-vector kn oc k± within the model. Numerical sim- 
ulations @, §, support the GS95 model. 

In this paper, we study compressible supersonic sub- 
Alfvenic MHD turbulence in low-/? plasmas. 



II. THEORETICAL CONSIDERATIONS 



While the GS95 model describes incompressible MHD 
turbulence well, no accepted theory exists for compress- 
ible MHD turbulence. Earlier theoretical and numerical 
efforts [[Io[ |2(], addressed effects of compressibility for 
limited parameter spaces. In (isothermal) plasmas, there 
are 3 types of MHD waves - Alfven, slow, and fast waves. 
Alfven modes are incompressible while slow and fast 
modes are compressible. Lithwick & Goldreich [jl6| con- 
jectured that Alfven modes follow the GS95 model and 
slow modes passively follow the same scalings for high (3 
(= P 9 /Pb ~ 2a 2 /V^; P g =gas pressure, PB=magnetic 
pressure; a=sound speed; V/i= Alfven speed) regime, 
which is largely similar to the exactly incompressible 
regime. They also mentioned that this relation can carry 
on for low (3 plasmas. 



In the ISM (3 is frequently less than unity. For instance, 
it is ~ 0.1 or less for molecular clouds. Therefore, we con- 
sider low [3 regime in this paper. Interstellar turbulence 
is traditionally thought to be sub-Alfvenic (SV < Va), 
although this is not a universally accepted assumption 
(see e.g. ||). If turbulence is super- Alfvenic initially, we 
expect that eventually magnetic energy should approach 
the equipartition level jl| and the scales smaller than the 
energy injection scale should fall in the sub-Alfvenic com- 
pressible regime. 



Arguments in GS95 are suggestive that the coupling 
of Alfven to fast and slow modes will be weak. Conse- 
quently, we expect that in this regime the Alfven cascade 
should follow the GS95 scaling. Moreover the slow modes 
are likely to evolve passively (see |li| ) , so that we expect 
the GS95 scaling for them as well. However, fast modes 
are expected to show isotropic distribution as their ve- 
locity does not depend on magnetic field direction. To 
test those theoretical conjectures we use numerical sim- 
ulations. 
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III. 



NUMERICAL METHOD 



To mitigate spurious oscillations near shocks, we 
combine two essentially non-oscillatory (ENO) schemes. 
When variables are sufficiently smooth, we use the 3rd- 
order Weighted ENO scheme 12 without characteristic 
mode decomposition. When opposite is true, we use the 
3rd-order Convex ENO scheme jl7| . We use a three-stage 
Runge-Kutta method for time integration. We solve the 
ideal MHD equations in a periodic box: 



dx/dt + v • Vv + p^V^p) 



dp/dt + V • (pv) = 0, 
(V x B) x B/471-p = f, 



dB/dt- V x (v x B) = 0, 

with V • B = and an isothermal equation of state. Here 
f is a random large-scale driving force, p is density, v is 
the velocity, and B is magnetic field. The rms velocity 5V 
is maintained to be approximately unity, so that v can 
be viewed as the velocity measured in units of the r.m.s. 
velocity of the system and B / ^Airp as the Alfven speed 
in the same units. The time t is in units of the large eddy 
turnover time (~ L/SV) and the length in units of L, the 
scale of the energy injection. The magnetic field consists 
of the uniform background field and a fluctuating field: 
B = B + b. 

For mode coupling studies (Fig. [I]), we use 144 3 grid 
points and we do not drive turbulence. We explicitly 
vary the Alfven speed of the background field, Va — 
Bo/y/ATrpQ, and/or the sound speed. Here p is the aver- 
age density. For scaling studies (Fig. |2|), we drive turbu- 
lence solenoidally in Fourier space and use 216 3 points, 
Va = 1, Po = 1) and a = VOX The average rms velocity 
in statistically stationary state is SV ~ 0.7. Therefore, 
the scaling results reported here utilize M s (= SV/a) ~ 
2.2, M A (= SV/Va) ~ 0.7, and (3 ~ 0.2. 



IV. RESULTS 

Mode coupling of MHD waves. — We first describe 
how to separate Alfven, slow, and fast modes in wave- 
vector (or, Fourier) space. In general, displacement vec- 
tors (hence v*,) of slow waves and fast waves are 



oc fc||k|| 



1 - VD - (3/2 
~JJ2 

-l 2 



D 



-i 2 



fcj_k_L, 



D + (3/2 



l + VD-p/2 



fc||k|| + fcj_kj 



(1) 
(2) 



where D = (1 + (i/2) 2 - 2(3 cos 2 9 and 8 is the angle 
between k and Bq. In the limit of (3 — > 0, the displace- 
ment vectors of the slow waves are almost parallel to k|| 
(|| Bo) and those of fast modes are almost parallel to 
k^ (_L Bo). We can obtain slow and fast velocity by 
projecting velocity component Vk into £ s and £/, respec- 
tively. We can obtain velocity and magnetic field due 



to Alfven modes in the same way as in incompressible 
case (see Q): £a = k|| x k^. To separate slow and 
fast magnetic modes, we assume the linearized continu- 
ity equation (tupk = Pok- v,t) and the induction equation 
(wbfc = kx (Bo x Vfe)) are statistically true. From these, 
we get Fourier components of density and non-Alfvenic 
magnetic field: 



Pk 



(p Av k ,s/c s )k ■ £ s + (p Av k j/cf)k ■ £f 
Pk,s+Pk,f, (3) 
(B Q Av Ks /c s )\B x | s | + (B Av k j/c f )\B x £/| 
h,s + b kJ (4) 

PkA B o/po){\B xii/k-i) 

p k j(B / Po )(\B xi f \/k-i f ), (5) 



where Av k oc 



(superscripts '+' and '-' represent 



opposite directions of wave propagation) and subscripts 
's' and T stand for 'slow' and 'fast' modes, respectively. 
From equations (||), (||), and (||), we can obtain pk, s , 
Pkji bk.s, and b k j in Fourier space. We obtain energy 
spectra (Fig. |^a,c,e) using this projection method done 
in Fourier space. When we calculate structure functions 
(Fig. |^b,f) wc first obtain the Fourier components using 
the projection and, then, we obtain the real space values 
by performing inverse Fourier transform of the projected 
components. However, we use a different method for the 
structure function of slow mode velocity (see Fig. |^d). 

The dispersion relation of Alfven modes and those of 
slow and fast modes in (3 — > limit are u) = Va^m, lu — 
ak\\, and to = Va&, respectively. Alfven modes are 
not susceptible to collisionless damping. Therefore, we 
mainly consider transfer of energy from Alfven modes to 
the compressible MHD ones (i.e. slow and fast). 

To check the strength of the coupling, we first per- 
form a forced supersonic sub-Alfvenic MHD simulation 
with Bo/y/inpo — 1. Using the same data cube obtained 
from this simulation, we perform several decaying MHD 
simulations. We go through the following procedures be- 
fore we let the turbulence decay. We first remove slow 
and fast modes in Fourier space and retain only Alfven 
modes. We also change the value of B preserving its 
original direction. We use the same constant initial den- 
sity po for all simulations. We assign a new constant 
initial gas pressure P g (26). Note that (3 — P g / (B / 8n) 2 . 
After doing all these procedures, we let the turbulence 
decay. We repeat the above procedures for different val- 
ues of Bq and P g . Fig. [l]a shows time evolution of kinetic 
energy of a simulation. The solid line represents the ki- 
netic energy of Alfven modes. It is clear that Alfven 
waves do not efficiently generate slow and fast modes. 
Therefore we expect that Alfven modes follow the same 
scaling relation as in incompressible case. Fig. [j]b shows 
that the following relation fits the data well: 



{5V) 2 f /{5V)\ ex (6V) A /B , 



(6) 



which means the coupling gets weaker as Bq increases 
p7|. Note that (SV) a and p are constants. This 
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FIG. 1: (a) Decay of Alfvenic turbulence. The generation 
of fast and slow waves is not efficient. Initially, f3 ~ 0.2 and 
Bo/y^Fp^ = 1. {b) The ratio of (SV) 2 f to (5V) 2 A . The ratio 
is measured at t ~ 3 for all simulations. The ratio strongly 
depends on Bo, but only weakly on (initial) /3. The initial 
Mach numbers span 1 — 4.5. 



marginal coupling agrees well with a claim in GS95, in- 
compressible simulations jl9|, and earlier studies where 
the velocity was decomposed into a compressible compo- 
nent and a solenoidal component ^ . 

Alfven Modes. — Fig. ||a shows that the spectra of 
Alfven waves follow a Kolmogorov spectrum: 



Alfven Waves: E A (k) oc fc ± 5/3 . 



(7) 



In Fig. |^b, we plot the second-order structure function 
for velocity (SF 2 (r) =< v(x + r) - v(x) > avg . over x ) 
obtained in local coordinate systems in which the parallel 
axis is aligned with the local mean field (see || ^, p^f). 
The SF2 along the axis perpendicular to the local mean 
magnetic field follows a scaling compatible with r 2 / 3 . The 
SF2 along the axis parallel to the local mean field follows 
steeper r 1 scaling. The results are compatible with the 
GS95 model (ry cx r 2 / 3 , or hi oc k 2 / 3 ). 

Slow waves. — The incompressible limit of slow waves 
is pseudo- Alfven waves. Goldreich & Sridhar || argued 
that the pseudo-Alfven waves are slaved to the shear- 
Alfven (i.e. ordinary Alfven) waves, which means that 
pseudo-Alfven modes do not cascade energy for them- 
selves (see also @ ) . We confirm that similar arguments 
are applicable to slow waves in low [3 plasmas. Energy 
spectra in Fig. ||c are consistent with: 



Slow Modes: E s (k) oc fc ± 5/3 . 



(8) 



In Fig. |^d, contours of equal second-order structure 
function (SF2), representing eddy shapes, show scale- 
dependent isotropy: smaller eddies are more elongated. 
The results are compatible with the GS95 model (fey oc 



fc 2 / 3 , or ry oc r 2 / 3 , where ry and r± are the semi- major 
axis and semi- minor axis of eddies, respectively ||). 

From the linearized continuity equation and the in- 
duction equation, we can show that density fluctuations 
are dominated by slow waves and only a small amount 
of magnetic field is produced by the slow waves in low 
fi plasmas: (Sp/p) s = (SV) s /a ~ M s , and (SB) S — ► 0, 
as (3 — ► 0. Here M s is the sonic Mach number. When 
M s ^> 1, the above relation for density fluctuation may 
not give a good approximation. 

Fast waves. — Fig. ||f shows fast modes are isotropic. 
The resonance conditions for the interacting fast waves 
are u)\ + u>2 = u>3 and ki + k2 = IC3 . Since u> oc k for 
the fast modes, the resonance conditions can be met only 
when all three k vectors are collinear. This means that 
the direction of energy cascade is radial in Fourier space. 
This is very similar to acoustic turbulence, turbulence 
caused by interacting sound waves Zakharov 
& Sagdeev [^4| found E{k) oc fc _3 ' 7 7 However, there is 
debate about the exact scaling of acoustic turbulence. 
Here we cautiously claim that our numerical results are 
compatible with the Zakharov & Sagdeev scaling: 



Fast Modes: E f (k) - fc~ 3/2 . 



(9) 



Non- Alfvenic magnetic field perturbations are mostly 
affected by fast modes when (3 is small: {8B)f ~ (SV)f, 
which is larger than (8B) S « 0. 

Turbulent cascade of fast modes is expected to be slow 
and in the absence of collisionless damping they are ex- 
pected to persist in turbulent media over longer times- 
pans than Alfven or slow modes. This effect is difficult to 
observe within numerical simulations where AB ~ Bq. 



V. CONCLUSION 

We found that, in the isothermal supersonic sub- 
Alfvenic low-/3 plasmas, the following scalings are valid: 

1. Alfven: E A {k) oc fc~ 5/3 , fey oc fc 2 / 3 , 

2. Slow: E s (k) oc fc~ 5/3 , fey oc fc 2 / 3 , 

3. Fast: E^(k) oc fc~ 3 / 2 , isotropic energy spectra. 
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